Geometrical and transport properties of sequential adsorption clusters 
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We investigate transport properties of percolating clusters generated by irreversible cooperative 
sequential adsorption (CSA) on square lattices with Arrhenius rates given by ki = q"' , where Ui is 
the number of occupied neighbors of the site i, and q a controlling parameter. Our results show a 
dependence of the prefactors on q and a strong finite size effect for small values of this parameter, 
both impacting the size of the backbone and the global conductance of the system. These results 
might be pertinent to practical applications in processes involving adsorption of particles. 



INTRODUCTION 

Transport properties in disordered media have been 
successfully and extensively studied within percolation 
theory [l|-|9| . Related to percolation is random sequential 
adsorption RSA of particles on a lattice [lO|. Sites are 
occupied sequentially independently of other sites and 
the system exhibits a spanning cluster at the percolation 
threshold. 

Despite the accomplishments of the RSA model in de- 
scribing disordered media [10, 11], interactions constitute 
an important feature in real physical systems. A common 
phase observed in adsorption is the c(2 x 2) |12l4l4|. in 
which particles are arranged in a checkerboard pattern 
[l5j , showing that nearest neighbor (NN) exclusion plays 
an important role in this process. This correlated filling 
has been studied for a chain molecule with reactive sites 
|16| and on square lattices, using both hierarchical rate 
equations 17[ and percolation theory [18'|, called in this 
case cooperative sequential adsorption (CSA). Filling of 
random and correlated 'lattice animals' has been the sub- 
ject of an extensive review [lOJ. The fractal dimension 
of the spanning cluster on CSA models has already been 
shown to be the same for traditional percolation |18l - 
|20[. However, the island formation process is known to 
be dependent on the filling ratios 2l|, |22| , and produces 
clusters with peculiar shapes for different ratios. Here 
we show and quantify the dependence of the transport 
properties of the spanning cluster on different ratios. 

In the CSA model, a particle irreversibly occupies an 
empty site i with a (multiplicative) rate given by 



h = g" 



(1) 



where rii is the number of occupied nearest neighbors of 
site i. The parameter q controls the strength an occupied 
site influences the filling of a neighboring unoccupied site: 
for g > 1 there is an attractive filling and for q < 1 a re- 
pulsive filling. For q — 1, every unoccupied site can be 
filled with the same probability, the model being equiva- 
lent to RSA. 
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FIG. 1. Representation of the CSA process on a square lattice 
for non interacting (g = 1), repulsive (g < 1) and attractive 
(g > 1) particles with the same occupied sites (in red). Empty 
sites are shown in grayscale. The probability of occupying an 
empty site increases with its lightness. For non interacting 
particles, the probability of occupying any empty site is the 
same, regardless the number of occupied nearest neighbors. 
For interacting particles it depends on the number of occu- 
pied NN. The difference between probabilities depends on the 
strength of the interaction, regulated by the parameter g, and 
the configuration of the system. 



METHOD 

We implement this model on a square lattice of size 
L with periodic boundary conditions in horizontal direc- 
tion. At each iteration step, we fill a random unoccupied 
site i chosen with probability 



P,:- 



q" 



Z-^Ui ^njQ' 



(2) 



where T„ . is the total number of sites with rij NN and 
q ^ 0- Fig. [T] illustrates the model for several values of q. 
For the special case oi q = 0, we notice that in the early 
stages of the process for very low q 



lim ToP(0) 



1. 



(3) 



Hence, only sites with unoccupied NN would be filled 
and the system would never percolate. Filling would oc- 
cur as in RSA with NN exclusion and stop at a iammed 
state with an occupation fraction of 0.365 17ll23l428l |. We 



modify the model for the special case of g — >■ 0+, labeled 
henceforth q = 0, to obtain results compatible with this 
limit in the following way. Sites are filled with NN exclu- 
sion until the system reaches the jammed state. Then all 
clusters in the system have size one, but no spanning clus- 
ter is formed yet. Next, sites with only one neighbor are 
then allowed to be occupied. When these sites are extin- 
guished, the ones with two neighbors are then occupied 
and so on, until a spanning cluster appears, promoting 
the global connection of the system. 

Once the spanning cluster is obtained, the backbone 
and the cutting bonds are identified using the burning 
method [29|. To calculate the conductivity, we assume 
that there is a resistor of conductance gtj between every 
pair of sites i and j of the system. The value of the 
conductance is 



9^3 



1, if i A j e backbone; 
0, otherwise. 



(4) 
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FIG. 3. Dependence on the parameter q of the prefactor in 
the finite-size scafing laws for the masses of spanning clus- 
ter, Ms = ObL'^'' , backbone. Mi, = a^L'^'' , and cutting bonds, 
Mc — acL^" ■ The cutting bonds prefactor Uc decreases mono- 
tonically with q but faster for attractive fiUing {q > 1). a, 
is close to ai, for repulsive filling, but for attractive filling it 
increases faster. 



An arbitrary current is then applied between two sites 
(i = 1 and i = Nback, the number of sites in the back- 
bone) in opposite sides of the spanning cluster and Kirch- 
hofF's current law is imposed at each site. As a result. 



we obtain a set of coupled linear algebraic equations: 



q = 



q = l 



q = 10 
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FIG. 2. Top: typical spanning clusters on the onset of per- 
colation for q = 0, g = 1 and q = 10. The spanning cluster 
appears in gray, the conducting backbone in blue and cut- 
ting bonds in red. For larger q, the fraction of sites of the 
spanning cluster present in the backbone is closer to 1. The 
number of cutting bonds generally decreases with q. Bottom: 
dependence of the percolation threshold Pc as a function of 
the parameter q. The orange line is the value for classical 
percolation, with g/(l + q) = 0.5. The dashed gray line is the 
value for q = Q. Average values are obtained for L — 4096 
and calculated over 2000 realizations. 






-1, ifi = l; 
1, iii = Ni,ack; (5) 

0, otherwise. 



where M is the Laplacian matrix of the backbone sites. 
Mij is given by 

Tii, if i = j; 
Mii = •^ — 1, if i and j are neighbors; (6) 

0, otherwise. 

We solve the system of Eqs. El for i,j = 2,3,..., N^us 
to obtain the voltages in each resistor and therefore the 
global conductance of the system. For each quantity an 
average is obtained over at least 2000 realizations for a 
given value of L. 



RESULTS 

The dependence on q of the percolation threshold [24| 
is calculated to high precision for system sizes up to 
L = 4096, as shown in Fig. [51 Also, typical realiza- 
tions of the system at the critical point are shown for 
three distinct values of q. Previous studies indicate that 
this model belongs to the same universality class of tra- 
ditional percolation [18|. Our data and analysis confirm 
that the fractal dimensions of the masses of the cluster, 
backbone, cutting bonds and of the conductivity are in- 
deed the same as in traditional percolation. 

Nevertheless, there is a clear difference in shap e for 
the clusters obtained for different values of q [18|. For 



very small values of q, the cluster resembles a warped 
wire and presents a large number of cutting bonds. As 
q increases and the repelling force becomes weaker (and 
turns to attraction for g > 1), the clusters become more 
massive. At the same time more sites are present in the 
backbone and there are less cutting bonds present. 

We measure this effect calculating the prefactors from 
power law fits at criticality for each mass with system 
size, Ms = asL'^% backbone, Mf, = a^L'^'' , and M^ — 
ttcL "", where dg, d^ and dc are the corresponding critical 
exponents. The results are shown in Fig. [3] While the 
spanning cluster mass prefactor, a^, varies little with q, 
the cutting bonds mass prefactor Uc decreases with q. 
For repulsive filling {q < 1) the backbone mass prefactor 
Of, is close to a„i but increases faster for attractive filling 
(g > 1). As a result, more sites of the spanning cluster 
are present in the backbone. These quantities, however, 
deviate from the expected power law behavior depending 
on the value of q. For g <C 1 there are strong finite size 
effects, as shown in Fig. ID We plot the masses of the 
spanning cluster, backbone and cutting bonds divided 
by the system size to the power of their corresponding 
fractal dimension (the prefactor for that system size) as 
function of q. For small g, the convergence towards the 
thermodynamic limit becomes very slow. We show that 
the low q prefactor converges to a finite limit Ooo as a 
power law of the system size. These values are shown in 
Table I. 

The results for the prefactor of the global electrical 
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FIG. 4. (a) — (c)Effect of the finite size of the samples on the 
scaling behavior of masses of the spanning cluster, backbone 
and cutting bonds. We show the mass divided by the system 
size to the power of the corresponding fractal dimension for 
different q. For small q, the convergence towards the thermo- 
dynamic limit becomes very slow, (d) Convergence of the floo 
towards the thermodynamic limit for q — >■ 0. We calculated 
aoo from the best power law fit of {M/L^f — a^c) for each 
quantity. 



TABLE I. Values of the low q prefactor a^c for which the 
curves in Fig. [3] converge, calculated by finding the best fit of 
{M/L''-f — Ooo) using a power law function. The exponent of 
the power law, a, is also given. 

Low q Prefactor aoo a 

Spanning cluster mass 0.4104 -0.99 ± 0.02 

Backbone mass 0.4708 -0.98 ± 0.02 

Cutting bonds mass 15.51 -0.51 ± 0.01 



conductance are shown in Fig. [S] for g G [0.1,4]. In 
this range it increases monotonically with g. For re- 
pulsive filling it increases slowly, being well fitted by a 
power law with exponent 0.13. For g = 4 the global 
conductance prefactor is more than 50% above the value 
of traditional percolation. This indicates that the more 
attractive the particles are, the better conductor is the 
conducting backbone cluster. 



CONCLUSIONS 

In summary, we have studied transport properties of 
percolating clusters generated by the CSA model with 
nearest neighbors correlations given by multiplicative 
rates modulated by a parameter g. Our results confirm 
that the mass of the backbone and the number of cutting 
bonds as well as the global conductance scale as power 
laws with the system size with the same exponents as 
in traditional percolation. We have shown that changes 
in the parameter g result in changes in all prefactors of 
these properties. These changes affect the spanning clus- 
ter shape, the coverage of the backbone, the number of 
cutting bonds and the global conductance of the system. 
As g becomes large, i.e., within the very attractive in- 
teraction regime, occupied sites are more likely to have 
occupied NN. This increases the fraction of sites in the 
backbone, allowing alternative paths for the current to 
flow in the cluster thus increasing the conductance of the 
system. Finally, we have shown the strong finite size ef- 
fect for g -C 1 and calculated the limiting values of the 
corresponding power-law prefactors. 
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